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Abstract 

Hydroxyl Proline dehydrogenase HYPDH is an enzyme used in the 
dehydrogenation of proline and is involved in the conversion of trans-4hydroxy 
proline to A-1- pyrroline-3-hydroxy-5-carboxylate oxidized to glutamate. 
PRODH2 is involved schizophrenia through microdeletion results in the loss 
of gene sequence. Consequently, the production of glutamate effected lead to 
a neurodegenerative disorder called schizophrenia which affects 1% of the 
population and is characterized by several behavioral abnormalities and lack of 
social behavior. Bioinformatics techniques were used for the prediction of 3D 
structures and protein-protein molecular docking studies of PRODH2. 
Threading and homology modeling approaches were employed for 3D 
structure prediction of PRODH2. STRING database was used and ALDH4A1 
was observed as the interacting partner of PRODH2. Protein-protein docking 
studies of PRODH2 against ALDH4A1 were done to analyze the interacting 
potential residues. Virtual screening was performed against FDA library from 
ZINC database and molecular docking was done by AutoDock Vina. It was 
observed that ZINC000016 molecule has least binding energy of -12.2 
kcal/mol. The scrutinized compound may have potential for further drug 
discovery processes. The observed results of this research may help in 
designing the novel therapeutic targets against PRODH2. 


© © This work is licensed under the Creative Commons Attribution Non- 
Commercial 4.0 International License. 
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Introduction 


Hydroxyproline dehydrogenase (HY PDH) is a protein 
encoded by proline dehydrogenase 2 (PRODH2) and 
also known as HSPOX 2. HYPDH belongs to the 
proline oxidase family [1, 2]. HYPDH performs 
numerous functions including use of flavin-adenine- 
dinucleotide to build a new electron transfer pathway 
and proline dehydrogenation. It involves the 
conversion of trans-4hydroxy proline to A-1- 
pyrroline-3-hydroxy-5-carboxylate oxidized to 
glutamate [3]. Glutamate is a wonderful 
neurotransmitter involved in signaling among neural 
cells [4, 5]. 

PRODH2 is known to express in many tissues 
however it has common expression in brain. It is a 
mitochondrial protein and involved in the transfer of 
redox potential. PRODH2 is present in the 
centromeric region of 22q11 and contains 14 exons. 
The 22qll region of chromosome undergoes 
microdeletion and known as a velocardiofacial 
syndrome (VCFS) [6, 7]. VCSF is characterized by 
distinct facial features and frequent infections. It has 
also been reported that the growth and development 
were also affected. 

The microdeletions occur at 22q11 leads to affect 
HYPDH results in a deficiency of the enzyme proline 
oxidase known as hyperprolinemic type 1 
characterized by high level of proline in the blood 
circulation [8]. Therefore, the production of glutamate 
reduced to perform the required function of 
neurotransmission as glutamate is the product of 
proline degradation pathway [9]. 

The deficiency of neurotransmitter results in various 
disorders as seizures and intellectual disabilities. It has 
reported that a high level of proline in blood 
circulation is associated with a low level of glutamate 
in the regions of brain that was weekly associated with 
schizophrenia a neurological disorder characterized 
by intellectual disabilities including hallucination 
[10]. 

Schizophrenia is a neurodegenerative disorder and is 
known to effect 1% of the population worldwide. It is 
a severe mental illness and also affects the family 
psychologically [11]. It is characterized by severe 
mental illness including delusions, disrupt speech, out 
of order thinking and lack of social behavior [12]. 
Schizophrenia is mainly caused by various genetic 
factors. Variations in many genes or deletions in many 
parts of the chromosomes results in the loss of 
information needed for a metabolism or pathway. The 
degradation of proline occurs and results in the low or 
completely lost expression of specific genes. 
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Additionally, with genetic factors there are many 
epigenetic and environmental aspects which lead 
towards schizophrenia [13]. 

It has reported that schizophrenia is associated with 
VCSF. The deletion of chromosome 22 undergoes 
microdeletions leads to genes lost. The results are 
severe and affect numerous metabolisms. The 
centromeric region of chromosome 22 encodes an 
enzyme called hydroxy] proline dehydrogenase takes 
part in proline degradation [11]. The microdeletion 
results in the loss of gene sequence and affects the 
production of glutamate results in weak 
neurotransmission and ultimately leads towards 
schizophrenia [14]. 

Computational drug designing have shown success in 
computational research methodologies to solve 
different biological problems [15] and 
computationally designed various molecules against 
cancer [16, 17] and neurological disorders [18-23]. 
The current effort has structural and virtual screening 
to identify novel compound against schizophrenia. 
The NMR and X-ray crystallographic structure of 
PRODH2 has not been reported yet. To predict the 
structural insight, a 3D structure of PRODH2 was 
predicted by utilizing computational 
approaches. 


numerous 


Materials and methods 


PRODH2 has one described isoform under the 
accession number of Q9UF12 in Uniport knowledge 
base and has six splice variants. In this study, 
sequence analyses, 3D structure prediction, virtual 
screening and comparative molecular docking studies 
have been performed. 

The amino acid sequence of PRODH2 was retrieved 
from Uniport KB (http://www.uniprot.org/) and 
subjected to blastp for the selection of suitable 
templates against Protein Data Bank (PDB) [24] for 
structure prediction. The MODELLER 9.20 [25] was 
employed to predict the three-dimensional (3D) 
structure of PRODH2. The comparative modeling is 
also performed by CPHmodels3.2 [26] which 
recognize templates based on profile-profile 
alignment and SWISS-MODEL [27]. Threading 
approach based programs including I-TASSER [28] 
and SPARKS-X [29] were also utilized. Some other 
structure prediction tools including Phyre2 [30], M4t 
[31] and Robetta [32] were also used to predict the 
structure of PRODH2. The 3D structure of PRODH2 
isoforms were visualized and minimized through 
UCSF Chimera 1.13 [33]. Anolea [34], ERRAT [35], 
Rampage [36] and Verify3D [37] evaluation tools 
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were used for the critical assessment of the structures 
and overall quality of the models. The Ramachandran 
plots were generated to evaluate ọ and y angles which 
showed the distribution of non-Glycine, non-proline 
residues. 

STRING (Search Tool for the Retrieval of Interacting 
Genes/Protein) database [38] was employed to 
analyze the interacting protein against PRODH2. The 
crystalline structure was obtained from PDB (PDB ID: 
4EO5). Gramm-X [39] online server was applied for 
protein-protein docking studies of PRODH2 isoform 
against 4EOS. PatchDock [40] was applied for further 
protein-protein interaction assessment. The 
hydrophobic and electrostatic interactions were 
mapped by using ligPlot [41]. 

Pyrx tool was used for the molecular docking analyses 
by optimizing AutoDock [42]. The interaction 
between the protein and ligand was employed for 
orientation and conformation. The FDA library of 
ZINC database was used for virtually screening of 
target protein against the library molecules. 





Results and Discussions 


Structural bioinformatics is one of the milestones 
archived by computational biologists. This dimension 
of biology aims to justify the macromolecular 
anatomy and physiology in a well-organized manner. 
The structural biology provides many disciplines in 
the field of pharmacoinformatics, 
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percussion medicine which helps to understand the 
potential of schizophrenia by analyzing various 
proteins or defective enzymes involved in 
schizophrenia susceptibility. Besides, all the 
researches holding in all over the World on exploring 
the effective remedies of schizophrenia and many 
researchers involved in it around the Globe. PRODH2 
is involved in schizophrenia as the dehydrogenation of 
proline by help in oxidization of L-proline to delta-1 
pyrroline-5-carboxylate. It was found that the 
variation in PRODH2 locus affects the availability of 
glutamate and mice with induced mutations exhibits 
abnormalities in their behaviors leads towards 
schizophrenia. In genome wide association studies, it 
was reported that the schizophrenia susceptibility 
locus was present at chromosome 22q11 region and 
microdeletion occurs leads towards schizophrenia. 


Sequence Analyses 


Extensive literature review and exploring biological 
databases revealed significant information about the 
sequence of PRODH2. UCSC genome browser was 
used to identify the exact location of PRODH2 in 
human genome. It was present on the q arm of 
chromosome 19 (Fig. 1). 

Clustal omega and NCBI blast [43] was used to get the 
MSA of the PRODH2 family to reveal the conserved 
regions for further analyses. Query coverage and 
percentage identity was analyzed having a word size 
of 3 and BLOSSOM62 matrix was utilized having a 
threshold of 10. 





Fig. 1: Location of human PRODH2 protein on chromosome 19. 


NCBI blast, ENSEMBLE and UCSC Genome 
Browser BLAT [44] were used for MSA. Clustal 
Omega was used for MSA of the sequences. Amino 
acids composition was calculated through Expesy 
ProtParam tool [45]. Numerous parameters were 
calculated including theoretical pI, which is 8.94, 
number of positively charged residues are 46 and 
negatively charged residues are 56 with the aliphatic 
index of 93.38 and the respective composition of 
amino acids were calculated (Fig. 2). 


Structure prediction 


Comparative modeling [46] and threading techniques 
were utilized for 3D structure prediction of PRODH2 
as the structure was not solved through X-ray 
crystallography and nuclear magnetic resonance 
(NMR). The sequence was retrieved in FASTA format 
from Uniport KB [47, 48] and subjected to blastp [49] 
to search for suitable templates against PDB for 
structure prediction. Three different templates having 
high query coverage, identity and E-value were 
selected as suitable templates for homology modeling 
having PDB ID SUR2, 4NM9 
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Fig. 2: The composition of the amino acid of PRODH2 calculated from ProtParam tool. 


and 5UX5 (Table 1). Modeller was used for 
homology modeling approach and models were 
generated against all the selected templates. It was 
observed that the maximum query coverage of the 
templates was 29% so threading approach it was 
utilized to lower down the error rate. 


Table 1: The best three templates 


Numerous structure predictions tolls including I- 
TASSER, SWISS-MODEL, M4T, phyre2, 
CPHmodels, SPARKS-X, Raptor-X and Robetta were 
used for implementing homology modeling and 








Description Max Score Query Cover __E Value Per. Ident Accession 
Chain A, Bifunctional protein PutA [Bdellovibrio 67.4 29% 3e-11 30.38% SUR2 
bacteriovorus HD100] 
Chain A, Proline Dehydrogenase And Delta-1-pyrroline- 52.0 24% 2e-06 29.23% 4NM9 
5-carboxylate Dehydrogenase [Geobacter sulfurreducens 
PCA] 
Chain A, BIFUNCTIONAL PROTEIN Proline utilization 43.1 29% 0.001 26.99% SUXS 


A (PutA) [Corynebacterium freiburgense] 





Note: Employed from the BLASTp 
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threading approaches against PRODH2. The predicted 
structures were further evaluated through various 
evaluation tools. Rampage tool was utilized to check 
the stereochemical parameters of the predicted 
structures. The phi and psi torsion angles were also 
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calculated. The favored region allowed region and 
outliers of the structures were calculated 

and the structures were further evaluated through 
ERRATE used to calculate the overall quality factor 
of the predicted models [56] (Fig. 4). 
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Fig. 3: Graph of quality factor, favored region, allowed region and outliers of the PRODH2 structures evaluated by 


using ERRAT and Rampage structure evaluations tools. 


All the predicted structures were evaluated on the 
values of their favored region, allowed region, outliers 
and overall quality factors. The comparative graph 
was generated from the evaluated values of the 
predicted structures. The most reliable structure was 
selected from the graphs for further analyses. 
Interestingly, it was observed that the highest overall 
quality factor was 100%. The minimization of the 


selected structure was performed by using UCSF 
Chimera 1.13 to improve the stereochemistry and 
optimization of the model through 1000 steepest and 
1000 conjugate steps. 


Protein-protein docking analyses 


String database was used to find the interacting 
partner of PRODH2, and it was observed that the 
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crystalline structure having PDB ID:40E5 was the 
interacting partner of PRODH2. The protein-protein 
docking studies were carried out by using GrammX 
online server and cross verify the results through 
PatchDock and FireDock. The interacting residues of 
the complexes were visualized and analyzed by using 
UCSF Chimera and ligplot. 

The interaction between PRODH2 and ALDH4A1 
(Table 2) were visualized and analyzed through 
UCSF Chimera (Fig. 6). ALDH4A1 was observed as 
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the interacting partner of PRODH2 leads to perform 
the protein-protein molecular docking analyses to 
analyze the residual interactions. The docked complex 
predicts the interacting. The protein-protein molecular 
docking studies were analyzed on the basis of least 
binding energy values. Numerous complexes were 
analyzed and top ranked complexes having least 
energy values for further refinement. The selected 
complex showed the least global binding energy and 
suggested that PRODH2 and ALDH4A1 have 
efficient binding affinity. 


Fig. 5: Protein-protein docking analyses of PRODH2 and ALDH4A1. 
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Compa rative molecular docking (Table 3). It was observed that the docked structure of 
PRODH showed reliable docking results. All the 


docked complexes showed similar binding domain. 
The top ranked ligand from selected library was 
elucidated, namely ZINC0000016 (Fig. 7) and the 
observed interacting residues were GLU504, 
ARGS500, ARU380, GLY381, ALA382, HIS436, 
VAL379, TYR497, LYS236, SER435, ALA434, 
GLN463, LYS377, LEU464, SER485, GLN234, 
ASP314 (Fig. 8). 


Virtual screening was performed for comparative 
molecular docking analyses of PRODH2 (Figure 7). 
The FDA library of ZINC database was screened by 
utilizing AutoDock Vina. Moreover, the FDA library 
query for virtual screening yielded 134 compound 
hits. The molecular docking analyses were done by 
AutoDock Vina. The top ranked least binding energy 
molecule was scrutinized for further analyses. The 
slight variation in least binding energy was observed 


Table 2: Interacting residues of the target protein PRODH2 and Interacting protein ALDH4A1. 








Target protein Target protein Residues Interacting protein Interacting protein residues 

PRODH2 GLY 103, LEU 373, ASP__ ALDH4A1 HIS53, SER50, ASP52, LEUS1, GLN307, HIS308, ASN155, 
52, VAL 215, TYR 162, VAL152, ASP102, GLY103, ARG104, LEU161, GLY339, 
GLY 452 TYR162, PROI63, LEU372, PRO427, GLY371, ALA423, 


GLY452, ILE453, HIS415, GLU39, ARG219, VAL215. 











A 


Fig. 8: Molecular docking analyses of PRODH2 and top ranked scrutinized ligand ZINC000016. 
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Table 3: Top ranked four compounds having least binding energy 








Compound name Binding affinity (kel/mol) RMSD (lower) RMSD (upper) 
ZINC000016 -12.2 0.0 0.0 
ZINC000059 -11.3 0.0 0.0 
ZINC000076 -10.8 0.0 0.0 
ZINC000060 -10.9 4.116 4.586 





Conclusion 


Schizophrenia is a complex disease that affects 1% of 
the population. In conclusion, in silico approaches 
were used based on structural modeling, virtual 
screening and molecular docking analyses. The 
predicted 3D structure suggested the reliability of the 
structure with good degree of accuracy. The molecular 
docking results satisfied the least binding energy and 
suggested that ZINC000016 may prove as a potent 
molecule for SZ by targeting PRODH2. The applied 
methodology leads to simplify the process of novel 
drug designing. 
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